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Abstract: We explain why naive discretization results that have appeared in [hep-lat/0006013] 
do not appear to yield the desired continuum limit. The fermion propagator on the lattice in- 
evitably yields a diagram with nonvanishing UV degree D = contribution in lattice perturba- 
tion theory, in contrast to what occurs in the continuum. This diagram gives a finite correction 
to the boson 2-point function that must be subtracted off in order to obtain the perturbation 
series of the continuum theory, in the limit where the lattice spacing a vanishes. Using a trans- 
fer matrix approach, we provide a nonperturbative proof that this counterterm suffices to yield 
the desired continuum limit. This analysis also allows us to improve the action to 0(a). We 
demonstrate by Monte Carlo simulation that the spectrum of the continuum theory is well- 
approximated at finite but small a, for both weak and strong coupling. We contrast the above 
situation for the naive lattice action to what occurs for the supersymmetric lattice action, which 
preserves a discrete version of half the supersymmetry. There, cancellations between D = 
diagrams occur, obviating the need for counterterms. 
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1. Introduction 

In this article we study lattice definitions of the simplest sort of Id supersymmetric quantum 
mechanics (SQM). 1 The continuum theory is invariant under two independent supersymmetries 
(susys) whose algebra is {Qi,^} = 2idt- From the imaginary-time continuum action, one can 
easily write down a naive lattice action. The next task is to study the quantum continuum limit. 
This exercise proves instructive as it illustrates certain issues and subtleties in the continuum 
limit of a lattice theory that contains fermions. The lessons that are learned have relevance to 
the lattice definition of supersymmetric field theories, which are the prime motivation of this 
work. 

The naive discretization breaks both of the supersymmetries of the continuum theory. A 
superior lattice action has been constructed by Catterall and Gregory (CG) [2]. Their action 

■"■For a review, and extensive references, see [1]. 
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possesses one exact lattice susy that is a discrete version of, say, Q\. However, for an interacting 
theory, this lattice action cannot be made simultaneously invariant with respect to a discrete 
version of Qi- The susy lattice action also follows from a general superfield analysis [3] of SQM 
actions that preserve 1 nilpotent supercharge, as well as the topological field theory approach 
of [4]. CG have studied the susy lattice theory by Monte Carlo simulation. Some results found 
by CG are of interest to us in this article. CG have computed the bosonic and fermion mass 
gaps. They find that the spectrum is degenerate, already at finite lattice spacing a. CG have 
compared these results to what occurs for the case of naive discretization. They find that the 
naive discretization gives a very different answer. Indeed, looking at their results, it would 
appear that the naive discretization does not yield the spectrum degeneracy of the continuum 
theory in the a — > limit. In particular, the bosonic mass gap does not approach the correct 
value in this limit. 

This result seems rather strange. CG have attributed it to large nonperturbative renormal- 
izations. In this article we take a closer look at the issue. We find that the modified fermion 
propagator that occurs on the lattice gives a nonvanishing diagram with UV degree 2 D = 0. In 
such a diagram, modes at the UV cutoff a -1 do not decouple in the a — > limit, and their effects 
need to be subtracted off with a counterterm. The susy lattice action is superior in this regard: 
an additional D = diagram appears, and just cancels the D = diagram of the naive theory. 
Consequently, UV modes do decouple in the a — > limit, and there is no need to introduce 
counterterms. These features explain the discrepant results of CG. 

We now summarize the remainder of this article. In Section 2, we review relevant aspects of 
the continuum theory. Then we describe lattice actions and the effect of discretization on susy. 
In Section 3, we discuss lattice perturbation theory and make our main points. The modifications 
to the UV behavior of perturbation theory that occur on the lattice have an important effect: 
one has to add a finite counterterm to ensure that a D = graph cancels. The lattice theory then 
has only graphs of negative UV degree; Reisz's theorem [5] guarantees that they approach their 
continuum value, leading to the desired supersymmetric continuum theory. The absence of these 
counterterms in [2] explains why CG did not obtain a susy continuum limit from the naive lattice 
action. In Section 3 we also present a transfer matrix analysis, which provides a nonperturbative 
proof that the counterterm suffices to guarantee the continuum limit. Furthermore, this analysis 
allows us to derive the C(a)-improved action. 

In Section 4, we describe Monte Carlo simulations of the naive lattice action, with and 
without corrections. We explore the effect of both (i) the necessary 0(1) subtraction and 
(ii) the 0(a) improvement that were derived in Section 3. We find that these corrections are 
sufficient to yield good agreement with the continuum at small but nonzero a. In Section 5 we 
draw our conclusions. 



2 Our definition of UV degree is the same as in [5]. An explicit description of UV degree will be given below. 
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2. Preliminaries 



In this section, we review elementary aspects of the continuum theory and set out our notations. 
Certain features of the naive lattice action and susy lattice action will be discussed. These 
considerations will set the stage for our later perturbative, transfer matrix, and Monte Carlo 
analyses of the quantum continuum limit. 

2.1 The continuum theory 

The hamiltonian of the theories that we study is of the form: 

H SQM = 1 -p 2 + 1 -h'\ q )- 1 -h"( q )[b\b} , 

[q,p]=h {b,b*} = l, 6 2 = (6+) 2 = . (2.1) 

Here, q and p are the bosonic coordinate and momentum operators, and b\ b create and destroy 
the fermionic state. That is to say, a basis of states is \x, ±) where 

q\x,±)=x\x,±), 6|a;,-) = 0, tf\x, -) = \x, +) . (2.2) 

We will refer to the function h(q) that appears in (2.1) as the superpotential; there, h'(q) = 
dh(q)/dq, etc. It is well-known that the theory is invariant under two independent susys whose 
algebra is {Qi,Q2} = 2H S qm = Hd t . These are nothing but Q\ = b(ip + h'(q)) and Q 2 = 
b\-ip + h\q)). 

It is straightforward to write a imaginary-time action corresponding to (2.1): 



dt 



o 



l -(x 2 + h' 2 (x))+^d t + h"(x))^ 



(2.3) 



Here, ip and ip are independent Grassmann fields, and x is real. All are one-component variables. 
The action is invariant under the imaginary-time continuation of the susy transformations, 
generated by infinitesmal Grassmann parameters 61,62: 

Sx — 6i*0 + 62*0 , 5$ — -ei(x + ti) , Sip = -e 2 (x - ti) . (2.4) 

To prove this invariance one need only make use of periodic boundary conditions for the fields. 
Equations of motion are not involved, which is important since they do not generally have a 
solution in the imaginary-time theory. With 5 = e\Q\ + 62(^2, 

Qix — ip, Qii> — -{x + ti) , Qiip = 0, 

Q2X = ij, Q 2 ^ = 0, Q 2 i/> = -(x - h') , (2.5) 

defines the supercharges Qi,Q 2 - 
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In what follows it will be convenient to distinguish as h the part of the superpotential that 
leads to interaction terms in the action: 

h = -mx 2 + h, h = Y —x n • (2-6) 

n>2 

A particularly simple case that we will concentrate on is the one studied by CG: 

h = ^mx 2 + ^gx 4 . (2.7) 

2.2 The naive action 

Using a naive discretization of the bosons, and Wilson fermions to avoid spectrum doublers, one 
obtains the following naive lattice action: 

„-S = \A- Xl A- x , + l KK + U^(r h + m . (2.8) 

(For convenience — here and below — we have moved the factor of a that comes from the dt of 
the continuum to the l.h.s.; repeated indices are summed.) Throughout this article, we denote 
the finite difference operators that we use according to: 

A s = -(A + + A"), A 2 = A-A+ = A+A" , (2.9) 

and A w (r) is the Wilson operator, defined in (2.11) below. We could have chosen the forward 
difference operator A + just as well in the bosonic part of the action; simple identities show that 
for periodic boundary conditions 

A'XiA'Xi = A + XiA + Xi . (2.10) 
As stated, we have introduced the Wilson operator 

A w (r) = A s - —A 2 (2.11) 

in the fermion kinetic term. In Id it interpolates between various choices for the finite difference 
approximation to the time derivative: 

A w {r = ±1) = A T , A^(r = 0) = A 5 . (2.12) 

Thus, the naive choices A + or A~ for the fermion kinetic term are implicitly contained in (2.8). 
Taking advantage of (2.10), we obtain the naive discretization of all time derivatives appearing 
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in (2.3) according to d t — > A - if r = 1 and d t — > A + if r = —1. These are free of doublers, 
whereas A s has doublers. There are two reasons why we use the more general operator A w (r). 

First, for d > 1, i.e. in the field theories that we ultimately have in mind, the naive choices 
of dfj, — > A^ (/i — 1, . . . , d) always lead to fermion doublers, as does A^. By contrast, the d- 
dimensional Wilson operator AYf(r) = A^ — (ra/2)A 2 does not have doublers, provided r ^ 0; 
it does not interpolate between the naive discretizations in d > 1. (Here, A 2 = A+A~ is the 
d-dimensional lattice laplacian.) Thus, we seek to understand the role of the Wilson operator 
in the simpler Id case as a foundation for d > 1 lattice susy studies. 

Second, the Wilson operator for r such that |r| <1 gives a "physical" interpretation to 
UV modes, as states of mass 0{ra' v ). This is particularly useful for our purposes in Section 3. 
We find a lack of decoupling of these "lifted doublers" in a diagram with UV degree D = 
(cf. Section 3.1). It is this effect that is responsible for a finite violation of susy Ward identities 
(cf. Section 3.3) — necessitating a finite counterterm to the naive action. Thus, we use the Wilson 
operator to illustrate the origin of the corrections coming from the UV scale. 

We choose to replace the continuum susy by a discrete approximation dt — > A + : 

QiXi = i>i, Qi4>i = ~{A + Xi + h'i), Q^i = 0, 

Q2Xi = A, Q 2 & = 0, Q 2 ipi = -(A + Xi - h'i) . (2.13) 

(Other choices could be made; they do not significantly alter our arguments. For example, 
it is possible to remove the quadratic terms in (2.14) below. However, it can be shown that 
the remaining terms will still lead to a finite violation of the susy Ward identities.) It is then 
straightforward to work out the variation of the action. For example, 

a^Q.S = ~(1 + r)^A-A 2 ^ - |(1 - r)ma:,A 2 ^ + y^A 2 ^ 

+ (A 5 ^ - hZA+XtWi . (2.14) 

In each equation we have substituted (2.6). Now we note that (no sum over i implied): 

A s ^ - h'fA+Xi = ~h'fA 2 Xi + 0(a 2 ) . (2.15) 

Thus we obtain that all terms on the r.h.s. of (2.14) are suppressed by a factor of a. (As is well 
known, this 0(a) violation is due to the failure of the Leibnitz rule on the lattice [6].) A similar 
statement holds for Q2S. The classical continuum limit therefore yields 

rP 

Q A S = ay A = aY] aY Ati -> I dt aY A (t) -> , A = 1, 2 . (2.16) 

Jo 

The susy invariance is recovered in this limit. On the other hand, in Section 3 we will find that 
in the quantum theory there is a finite violation of the continuum susy Ward identities: 

(Q A S) = a(y A ) ~ a ■ a' 1 -> finite . (2.17) 

The main point of this article is to explain why this occurs and how to cure it. 
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2.3 The susy lattice action 

Here the action is the one that appeared first in [2]. It preserves a discrete version of half the 
continuum susy exactly. The action was subsequently rederived in [3] using a lattice superfield 
approach, and in [4] using a topological field theory approach. It is given by (we suppress the 
r-dependence of A w ): 

a-'S = \ £(A^ + h'f • + . (2.18) 

i ij 

Comparing to the naive action (2.8), it can be seen that the fermionic part of the action is 
unchanged. By contrast, the bosonic part of the action differs from the naive discretization (2.8) 
in a number of ways. Due to supersymmetrization, the Wilson operator A w appears in both the 
boson and fermion kinetic terms. It follows that at r = the bosons also have doublers; these 
are similarly lifted by a Wilson mass term m w = m — (ra/2)A 2 for r ^ 0. This degeneracy of 
UV modes is necessitated by the one exact lattice susy. The lattice susy also requires derivative 
interaction terms that are not present in the continuum theory. To see this, note that the 
bosonic part of the action Sb may be written as: 

a- l S B = i^[-A 5 A s + (m - yA 2 ) 2 ]^ + mtipi + ^hfa 

+h' i A s x i - ^h'^Xi . (2.19) 

It is clear from (2.19) that in addition to the modified quadratic action, there are also (in 
the second line) two sorts of interaction vertices not present in the continuum. The vertices 
associated with h' i A s x i come from the h! i A w Xi crossterm one gets from the square in (2.18). 
As pointed out in [3], these interactions break reflection positivity; nevertheless the continuum 
hamiltonian and Hilbert space are recovered after a "conjugation" of states and the transfer 
matrix — by the operators exp(±h(q)), where q\Xi) = Xi\xi). It should also be noted that whereas 
f h'dx = f dh vanishes in the continuum with periodic boundary conditions, on the lattice 
h' i A s Xi is not purely a total derivative, due to the failure of the Leibnitz rule. The vertices 
associated with rah' i A 2 Xi are "Wilson interaction terms." Both of these sets of vertices have 
UV degree D — 1, so that lattice power counting differs from that of the continuum. We will 
examine this issue in Section 3.4 below. 

As mentioned above, our detailed considerations will focus on the superpotential (2.7). In 
what follows, we find it convenient to set r = 1; this leads to various simplifications of the above 
formulae, which we summarize here: 

M ik = a(A~ + h" k ) = -5i_i,fe + fii6 ik , 

Hi = l + a(m + 3^ 2 ), detM = -1 + JJ ^ . (2.20) 
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The matrix M is the one which appears in the fermionic part of both the naive and susy actions: 
Sf — ipiMikip- For m > and g > 0, the fermion determinant is strictly positive. This is a 
happy state of affairs, given the sign/complex-phase problems that are often experienced in susy 
lattice systems; for example [7]. We note, however, that when h is a polynomial of odd degree, 
positivity will not hold. (This happens to be the case where it is known from the continuum 
that susy is spontaneously broken.) 

As first pointed out by CG, the action may be written in the Nicolai map form: 

a^S = ^MM + , Mi = A w x t + h'i . (2.21) 

This form makes the exact susy rather obvious: 

5 x Xi = e^i , S^i = , S^i = -eiMi ■ (2.22) 

The index "1" indicates the correspondence to the e\ part of the continuum susy (2.5), the 
lattice version of Q\ which has been preserved exactly on the lattice. It is easy check that this is 
also a symmetry of the partition function measure. This comes with the important qualification 
that periodic boundary conditions are always assumed for both bosons and fermions. Thus 
our partition function is Z = Tr (— l) F e~ l3H , the Witten index. Continuum (imaginary-time) 
vacuum expection values are obtained in the limit (3 — > oo. 

3. Counterterm analysis 

We will show below (Section 3.1) that the naive lattice action does not give rise to a finite 
lattice perturbation series, in the technical sense that a UV degree D = proper vertex exists. 
This is in contrast to the continuum (D — contributions vanish) and the susy lattice (D — 
contributions cancel — see Section 3.4). A local counterterm can be introduced to render the 
naive lattice theory finite. 3 We determine this counterterm and show by various means that the 
subtracted theory has the correct continuum limit. In the process, we are able to compute the 
C(a)-improved action (cf. Section 3.2), which has a faster convergence toward the continuum 
limit. In Section 3.3, we illustrate how the continuum Ward identities are recovered in the a — > 
limit, as a consequence of the finite counterterm. 

The D = diagrams that occur are just 1-loop corrections to the boson propagator. D > 
primitive diagrams do not occur at higher loops. We will find that the a — > results in the naive 
lattice case differ from that of the continuum due to modes at the edges of the Brillouin zone. 
Thus the discrepancy can be associated with fermion doublers that are lifted by the Wilson 
mass — the contribution of the second term in (2.11) at the edges of the Brillouin zone. This 

3 The spirit of this analysis is quite similar to that of [8], where a super-renomalizable supersymmetric field 
theory was obtained from a lattice theory by the treatment of counterterms. 
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non-negligible effect of doublers should not come as a surprise: the effects of lifted doublers, 
with 0(a~ l ) mass, need not decouple in a D = diagram. Rather, we find that they contribute 
a finite result. (While afl = diagram could contribute a term ~ In a, we will find that this is 
not the case and only an a- independent counterterm will be required.) 

As will be shown in Section 3.4, in the susy lattice action, an additional D = diagram 
appears. It cancels the finite correction that appears in the naive lattice theory, alleviating 
the need for a counterterm. The additional diagram arises because of the boson interactions 
associated with supersymmetrization of the Wilson mass term. Thus one advantage of the susy 
lattice action is that it provides a lattice perturbation theory that is finite (in the technical 
sense — only D < contributions survive for all proper vertices). By the theorem of Reisz [5], 
the susy lattice action therefore requires no counterterms to achieve the desired continuum limit. 
This explains why it was previously found [2] that the susy lattice action extrapolated well to 
the continuum, whereas the naive lattice action did not. 

We then show in Section 3.2 that the counterterm that is required for the naive action can 
be seen very easily in the transfer matrix description utilized in [3]. Remarkably, we are able 
to determine the counterterm that is required for an arbitrary superpotential, without having 
to perform any loop calculations. With a little more algebra, we are able to obtain the 0(a)- 
improved action. The transfer matrix analysis has the advantage that it is nonperturbative. 
Thus we are able to conclude that the counterterm suffices to guarantee the continuum limit be- 
yond perturbation theory. This is consistent with our simulation results, which will be discussed 
in Section 4. 

3.1 Power counting for the naive lattice theory 

As elsewhere, our focus here will be on the case of the superpotential (2.7). We will include 
all quadratic terms, including the Wilson mass for the fermion, in the free propagator. This 
does not alter the UV degree of internal fermion lines: since (ra/2)A 2 ~ a -1 at the edges of the 
Brillouin zone, the Wilson fermion propagator has the same degree 4 D = 1 as A s ~ a -1 . The 
interaction vertices can be worked out from (2.8): 

a-^f • = (™9*\ + \fA + 3^¥^) • (3.1) 

We count the number of x 4 vertices by V4, the number of x 6 vertices by V 6 , and the number of 
x 2 tptp vertices by V4. Each of these ver tices has UV degree D — 0. Then in the usual manner a 
4 The lattice degree of divergence of a lattice Feynman diagram is defined as the exponent D in the scaling 
a~ D of a lattice Feynman graph at small lattice spacing. In the present case of d — 1, a factor a -1 is assigned to 
every loop momentum integral (finite sum). The small-a scaling of bosonic (a 2 ) and fermionic (a 1 ) propagators 
is read off the expressions for lattice propagators. The scaling of vertices is read off from the vertex function 
that appears in the Feynman rule. For example, none of Feynman rules for the interactions in (3.1) scale with 
a; thus they have degree D = 0. On the other hand, the two new interaction in (3.22) below have Feyman rules 
~ ga^ 1 sin(2nk/N) ~ a -1 and ~ rag ■ a~ 2 sir? {irk / N) ~ a -1 resp. Thus they have degree D = 1. 
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diagram will satisfy the following relations between the UV degree D, the number of loops L, 
the number of internal boson and fermion lines I B , If resp., the number of external boson and 
fermion lines E B ,E F resp., and the various vertices: 

D = L - 21 B - I F , L = 1 + I B + I F - V A - V 4 - V 6 , 
E B + 2I B = 4V 4 + 6V 6 + 2V A , E F + 2I F = 2V 4 . (3.2) 

The only D > solution with L > is: 

D = E F = I B = V A = V 6 = 0, L = I F = V 4 = 1, E B = 2 , (3.3) 

the 1-fermion-loop correction to the boson 2-point function, diagram (a) in Figure 1. In the 
continuum, the corresponding expression is: 



dp —ip + m ( 1 



s - = 69 Zt Stt^ = 69 It + 0(ma) ) ■ {3A) 

(Note that we take the UV cutoff to be A = 7ra -1 , as is usual when relating lattice perturbation 
theory to the continuum.) We notice that the D — part of the diagram (i.e., involving the 
— ip part of the numerator) vanishes because it is odd w.r.t. p — > —p. For this reason, the result 
is finite in the a — > limit, rather than log-divergent as we would naively expect based on 
power-counting alone. At the same order in perturbation theory, a 1-scalar-loop diagram (b) of 
Fig. 1 must be added to this; it differs only by a factor of —2 relative to (3.4), yielding the net 
result -3$(1 + 0(ma)). 

'':/' 




(a) (b) 



Figure 1: Diagrams contributing at O(g) to the boson propagator. Diagram (a) has only a D = — 1 
piece surviving in the continuum, but on the lattice has a D = contribution coming from fermion 
doublers. In the susy lattice action, additional interactions due to supersymmetrization of the Wil- 
son mass term cause diagram (b) to also acquire a D = contribution, which just cancels that of 
diagram (a). 

On the lattice, this continuum expression is replaced by: 

6g -ia' 1 sin(2irk/N) + m + 2ra~ 1 sm 2 (7rk/N) 
latt - ~ iV^ a~ 2 sm 2 (2irk/N) + (m + 2m- 1 sin 2 (vrA;/A^)) 2 ' (3 ' 5) 
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A naive a — > limit of this expression would (i) set p = 2nk/(aN), (ii) replace (1/Na) J2k ~^ 
f (dp/2ir), and (iii) take the a — > limit, with the physical momentum p = 2irk/(aN) fixed, on 
the summand that appears in (3.5). Then one of course recovers the continuum fermion loop 
result (3.4). However, if the a — > limit is taken after computing the sum, one finds that the 
fermion loop contribution to the boson 2-point function is twice as big as (3.4): 

^latt. =6g(l + 0(rma)) . (3.6) 

This result can be understood as follows. The — ia~ x sm.{2nk/N) ~ — ip term is odd about 
k — and vanishes under the sum, just as in the continuum. However, the Wilson mass 
term 2ra~ l sin 2 (irk/N) is even and does not. Using p = 2nk/(aN), it can be written as 
2ra~ 1 sin 2 (pa/2). In step (iii) of the naive limit described above, we would neglect this term. 
However, in the explicit evaluation of the sum (3.5), the Wilson mass term becomes 0(a^ 1 ) at 
the edges of the Brillouin zone, where k ~ N/2. In order to take this non-negligible contribution 
into account, step (iii) should be avoided, 5 and the a — > limit must be taken after computing 
the sum. Since the Wilson mass term scales like a -1 , the naive expectations based on power 
counting are correct in the lattice theory — unlike in the continuum theory — there is a regulator 
dependent contribution that must be subtracted off. (However, a log-divergent result is not 
obtained for this particular D = diagram.) Whereas in the continuum we obtain (3.4), on the 
lattice there is an extra (finite) contribution of 1/2 coming from the edges of the Brillouin zone. 
The doublers do not decouple, precisely because this is a D = diagram. 

In order to match the continuum theory, we must add a boson mass counterterm that 
compensates for the effects of the doublers in the above diagram: 

_% v- 1 2ra-'sm 2 (irk/N) ^ 

Na^ o a~ 2 sm 2 (27ck/N) + (m + 2ra- 1 sm 2 (iik/N)) 2 ^ 9 m ■ ■ 
The shifted action is then 

s < = s + \ E W + const -] = s + \ E h " > ( 3 - 8 ) 

i i 

where we have added an overall constant (m/2) to obtain the appealing form of the far r.h.s. 
It will turn out that this constant appears automatically in our transfer matrix analysis of 
Section 3.2 below. We have only shown that a mass counterterm is required for the special case 
of h studied here, eq. (2.7). The expression (3.8) would, for a more general h, imply interaction 
counterterms as well. In Section 3.2 we will show that the counterterm h"/2 is the correct choice 
in the general case. 

We now give an intuitive understanding of why a finite counterterm is needed. Our dis- 
cussion will also give the more careful treatment of the region where p = 0(a~ 1 ), mentioned 

5 Or, as we show below, one can be more careful with the region where p = 0(a^ 1 ). 
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in Footnote 5, and will exploit the flexibility allowed by our use of the Wilson operator (2.11). 
Consider the case of r and a such that r < 1 but m <C ra^ 1 . Then the sum (3.5) is well 
approximated as follows. The mass of the (lifted) doubler mode is just M = 2ra~ 1 + m, as can 
be seen by a change of variables to p' = p — na" 1 , and looking near small p'. Thus for M <C a^ 1 , 
the doubler contributes like a particle of mass M ^> m. It follows that in this limit a good 
approximation to the integral that appears in (3.6) is just 



n/a dp m f w/a dp' M 1 r 



w/a 27rp 2 + m 2 J_ n/a 2tt P ' 2 + M 2 tt 



_i 7T , 7T 

tan h tan 



2am 2aM 



(3.9) 



r/a ^" V ' " b J-7r/a 

To 0(r), each term on the r.h.s. of (3.9) yields 1/2. Thus as the UV cutoff a -1 is sent to infinity, 
the contribution of the doublers does not decouple. For D > diagrams this lack of decoupling 
is usually associated with a divergence. Here, the lack of decoupling in a D = diagram is 
instead only associated with a finite, a-independent contribution coming from the regulator. 

Finally, we note that our result is consistent with the power-counting theorem of Reisz [5]. 
In essence, his theorem states that the naive continuum limit can be taken on the integrand of 
diagrams — step (iii) of the naive limit mentioned above — provided the UV degree of the diagram 
satisfies D < 0. That is, for D < diagrams, the a — > limit and the loop momenta integrals 
commute. However, this theorem does not apply to a D = diagram, regardless of whether or 
not it is finite (in the sense of being a-independent versus In a dependent). 

3.2 Transfer matrix analysis of the naive and improved lattice actions 

The partition function Z na i ve of the naive lattice action (2.8), with r = 1, can be represented as 
the trace of the iV-th power of the transfer matrix operator, denoted here by T[0,0], over the 
Hilbert space (2.2), weighted by (— 1) F to account for periodic boundary conditions of fermions. 6 
Thus, as we show in some detail in the appendix: 

Z nawe = Tr(-l) F T[0,0} N , (3.10) 

where T[0, 0] = T[k = 0, £ = 0] can be read off eqn. (3.12) below. In anticipation of counterterms, 
we generalize the naive action by functions k(x) and £(x): 

h' 2 {x) -> h' 2 {x) + k(x) , h"(x)^ -> (h"(x) + £(x))4>4> , (3.11) 

and denote the corresponding "deformed" transfer matrix as T[k,£], given by the following 
operator: 

oo 

T[k,£]=Af(a) J dzexp - ^h' 2 (q) - ±k(q) e iz *> [l + a(h"(q) + £{q))tfb] , (3.12) 



6 The case of r ^ 1 can also be treated, but it is more complicated because the transfer matrix must have row 
and column indices labeled by pairs of sites. 
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where J\f(a) is an inessential a-dependent normalization constant. (In the appendix we show 
that the naive lattice partition function supplemented by the counterterms (3.11) can be written 
as Tr(— 1) F T[k, £] N .) 

For small lattice spacing, one can use a saddle point approximation to evaluate the z-integral 
in (3.12). Then, it follows that the operator T[k,£] of (3.12) can, for small a, be written as: 

T[k,£] = e afiq) e - aH[k ' e] e- afiq) , (3.13) 

where / is given by: 

/(?) = -\(h' 2 {q) + k(q)) - \{h"{q) + £(q)pb , (3.14) 
and the hamiltonian H[k,£] is: 

H[k,£] = \ V 2 + 1%)' 2 - (h"(q)+£(q)) (l - 1 -a(h"(q)+£(q)^ \[b\b] 

+ 1 -k(q) - l -{h"{q) + £(q)) (l - l a {h\q) + €(</))) + 0(a 2 ) . (3.15) 

Two important observations immediately follow from (3.13-3.15): 

1. The T[k,£] operator at small lattice spacing is equivalent, upon conjugation by e a f , to 
exp(— aH[k, £]). When inserted in the trace, the factors e ±a ^ in (3.13) do not affect the 
partition function. (Actually, this conjugation can be avoided if one writes T[k,£] in the 
standard form e~t y ('?,fe t fe) e -fp 2 e -f v(q,tfb)^ w ith an appropriate choice of V .) Thus, the 
general partition function Z[k,£] = Tr(-1) F T[A;, £} N approaches Tr(-) F e^ H ^ in the 
small-a, /5-fixed limit, where (5 = Na. 

2. The hamiltonian H[0, 0] — corresponding to the naive lattice action (2.8) — does not, even 
as a — > 0, approach the hamiltonian Hsqm (2.1) of SQM. Rather, as eq. (3.15) shows, 
H[0,0]=H SQM -lh"(q) + O(a). 

The second observation above shows that in order for the continuum limit of the naive lattice 
partition function to correspond to the path integral of SQM, the naive lattice action has to be 
supplemented by a finite counterterm: one must take k(q) = h"{q) and l(q) = in the transfer 
matrix (up to 0(a) terms). This counterterm, with h of eq. (2.7), is precisely equal to the 
counterterm calculated in perturbation theory in the previous section, eq. (3.8). 

Moreover, eq. (3.15) allows us to go further and demand that the correspondence to H S qm 
is not corrected also at 0(a). One chooses k and £ such that H[k,£] = H S qm + 0(a 2 ), which 
implies 

k - (h" + £) ^1 - \{b!' + £)^j = 0(a 2 ) , 

h" - (h" + £) (l - ^a(h" + £)) = 0(a 2 ) . (3.16) 
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Thus, the 0(a)-improved naive lattice action should have a transfer matrix T[k,£] with 

k(q) = h"{q) + 0(a 2 ) , £(q) = \ah" 2 {q) + 0{a 2 ) . (3.17) 
The (9(a)-improved action is therefore (recall we set r = 1): 

a ~ ls ca = \Y. ( a ~^ a ~^ + KK + K) + ^(A". + h'fSij + \ati?ki)i>i ■ (3-18) 

i ij 

The above analysis provides a complete, nonperturbative proof that the 0(\) correction 
(3.8) to the bosonic part of the action suffices to guarantee the desired continuum limit. The 
O(a) improvement, which only corrects the fermionic part of the action, will be shown in our 
simulation results of Section 4 to yield very good results in the a — > limit. 

3.3 Susy Ward identities in the naive lattice theory 

It is easy to show that the following Ward identities hold on the lattice: 

(QaO) = ((QaS)O) = a(y A 0) , A =1,2. (3.19) 

Here 3^4 are the lattice operators defined by (2.16), whose explicit form can be obtained from 
(2.14) and a similar expression for Q 2 S. A violation of the continuum Ward identity (QiO) = 
can only arise if the difference operators in y 1 generate a UV-divergent result, (3^0) = 
0(a~ n ), n > 0. We now show that this indeed occurs in the uncorrected theory. This is in spite 
of the fact that the lattice perturbation series has only af = diagram, as shown in Section 
3.1. The point is that, as can be seen from (2.14), has vertices that have positive UV degree, 
whereas those that occur in the action (3.1) have vanishing UV degree. We will also illustrate 
how this difficulty is cured by the counterterm (3.8). 

We note from (2.14) that each term in contains finite difference operators. When these 
operators act on would-be fermion doubler modes (i.e., doublers when r = 0), they give positive 
UV degrees. We will study the anomalous Ward identity associated with the operator x^j, 
since it is directly related to the lack of spectrum degeneracy in the unsubtracted naive case. 
Actually, we must be careful to restrict to the physical modes that are relevant to the "long- 
distance" physics. Operators at fixed time could easily be obtained by blocking the lattice fields, 
but it is simpler to just analyze the Ward identies in momentum space; i.e., using the Fourier 
transformed variables Xk,ipk- Thus we will consider the Ward identities associated with the 
operator Ok = Xk^Pk subject to k N/2 mod N. We will find that whereas (Q\Ok) does not 
vanish in the continuum limit for the naive action, once the counterterm in introduced, the finite 
violation is subtracted off and the Ward identity is recovered in the a — > limit. 

Rather than study this cancellation directly, we will exploit the r.h.s. of the lattice Ward 
identity (3.19). We break Q X S into operators that generate 2-point and 4-point vertices: 



QiS = (QiS)(2) + (QiS)(4) = ay lt (2) + 0^1,(4) 



(3.20) 



The operator 3^1,(2) gives a D = 3 vertex and 3^1,(4) gives a D = 2 vertex. The leading order 
diagram associated with the latter vertex is shown in Fig. 2(a). The leading D = 2 piece 
vanishes. The surviving D = 1 part yields an (9(a _1 ) result that just cancels the factor of a 
in (3.20), leading to a finite result of 3g(l + 0(ma)). The 2-point vertex also contributes at 
O(g), through the diagram of Fig. 2(c). Here, the 3gx 2 if} , ifj vertex of the action (3.1) is involved. 
The diagram has UV degree D — 1, again canceling the a in (3.20), yielding a finite result 
of —6(7(1 + 0(ma)). The net result is finite and nonzero, which is consistent with the lack of 
spectrum degeneracy observed in [2]. The continuum Ward identity (QiOk) = is not recovered 
in the a — > limit, for the uncorrected naive action (2.8). 

However, once the counterterm (3/2)gx 2 is added to the action, 

QtS ^ Qrf + a^ZgxuPi . (3.21) 

i 

This contributes a 2-point D = vertex that is not suppressed by a, and is indicated by the 
tree-level diagram in Fig. 2(b). This of course yields 3g, canceling the finite violation that arises 
from the other two diagrams. Thus in the a — > limit, we obtain to O(g) the continuum result 
(QiOk) = 0. We have verified that similar cancellations occur at 0(g 2 ). 

Of course, it is guaranteed that the cancellations will occur to all orders in perturbation 
theory, since the counterterm has ensured that we will recover the continuum perturbation se- 
ries, which satisfies the susy Ward identities to all orders — for any operator O, not just the one 
that we have considered here. Furthermore, our transfer matrix analysis provides a nonpertur- 
bative proof that the continuum hamiltonian is obtained in the quantum continuum limit. We 
necessarily recover the continuum result (QaO) = in the corrected lattice theory when a — > 0. 




(a) (b) (c) 

Figure 2: Diagrams associated with the cancellation of the O(g) violation of the continuum susy Ward 
identity. The 2-point and 4-point shaded vertices that violate fermion number are those arising from 
QiS c , where S c is the corrected action indicated by (3.8). The sum of these diagrams vanishes in the 
a — > limit, provided the external momentum satisfies \p ex t\ < a" 1 . 

3.4 Power counting for the susy lattice theory 

In the transfer matrix analysis of [3], it was already proven for the susy lattice action (2.18) 
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that the correct continuum limit is obtained without the need for counterterms. However, it is 
interesting to see the detailed cancellations that bring this about in perturbation theory. 
In the special case (2.7), the interaction part of the action is given by: 

a-'SZT = E ( m 9< + \g 2 < + gx^A s Xi - "-fx^x, + 3gx^?j . (3.22) 

Two new terms appear in (3.22), relative to the naive lattice theory (3.1). Both of these vertices 
scale like a -1 and thus contribute to D. We include these into the definition of the 4-point boson 
vertex, associated with V 4 in our power counting analysis. We then obtain a slightly modified 
version of (3.2): 

D = L + V 4 - 2I B - I F , L = 1 + I B + I F - V 4 - V 4 - V 6 , 

E B + 21 B = AV A + 6V 6 + 2V A , E F + 21 F = 2V 4 . (3.23) 

It is straightforward to show that the only solutions with D > 0, L > are just (3.3) and 

D = E F = I F = V 4 = V 6 = , L = I B = V 4 = 1 , E B = 2 , (3.24) 

corresponding to the scalar loop Fig. 1(b). Thus in the susy theory we have two D = 1-loop 
diagrams to sum in the proper vertex associated with the boson 2-point function. It suffices to 
check at vanishing external momentum p ext , as higher orders in the expansion about p ext = 
will have D < 0. One obtains from the boson loop diagram 

v- 1 2m + 2ra- 1 sin 2 (nk/N) 

Na a' 2 sm 2 (2nk/N) + (m + 2ra~ 1 sin 2 (7rfc/iV)) 2 ' 1 ' ' 

Added to the fermion loop diagram (3.5), the D = contributions cancel. This is a consequence 
of the 1 exact susy of the lattice action. What is left over is a D = — 1 result. By the theorem 
of Reisz [5], as a — > this is guaranteed to give the continuum result — a finite mass correction 
—3g. (In the present case, it is trivial to check this explicitly.) 

4. Simulation results 

All of our simulation results correspond to the special case of superpotential (2.7) and r — 1. 
We will take m = 10, for the following reason. The effects of interactions always cause the mass 
gap mi (energy of the first excitation) to satisfy mi > m. Thus the Compton wavelength of this 
mode is always shorter than 1/10. For this reason we can safely set the system size f3 = Na to 
(3=1: the finite size effects should be negligible since several Compton wavelengths fit within 
the total volume, giving a good approximation to the (5 — > 00 limit. 
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All of our simulations are performed using hybrid Monte Carlo techniques. See [2] for 
details relating to the implementation for lattice SQM systems. Unlike [2], we do not use 
Fourier acceleration; rather, we have measured autocorrelation times and simply spaced our 
samples appropriately. Due to the simplicity of the fermion matrix, we were able to invert 
M T Ms = <j) analytically, where is a pseudofermion configuration, using two steps of Gaussian 
elimination. Alternatively, since M T M is a cyclic tridiagonal matrix, we were able to apply the 
Sherman- Morrison algorithms with partial pivoting. In either case, for the double precision that 
we used, we were able to keep the relative residual error \(M T M) ■ s est — (/>\/\<f>\ of the inversion 
to less than 1 x 1CT 14 throughout the simulations. This is far more accurate than what was 
obtained using conjugate gradient. Finally, fit errors due to excited states were avoided by only 
fitting times t of Green functions G(t) such that these effects could be shown to be negligible. 

We have extracted excitation energies, or, effective masses, from connected Green func- 
tions G J (t) where t = a, 2a, . . . , Na is the imaginary-time of points on the lattice, and / G 
{IB, IF, 2B, 2F} labels which excitation dominates the large-time behavior of the Green func- 
tion. In particular, we choose 

G lB (t) = (x lXl+t/a ) , G 1F (t) = (M 1+t/a ) , 

G 2B (t) = {xlx\ +t/a ) con ^ , G 2F (t) = (xnkXl+t/atl+t/a) ■ (4-1) 

Due to the symmetry x — > —x of the action, as well as fermion number, the states that contribute 
to each of these Green functions come from different sectors of the state space. For t <C Na and 
N >> 1, we have for example 

G lB {t) = c lB e~ m ' Bt + c 3B e~ msBt + ... , G 2B (t) = c 2B e- m2Bt + c 4B e~ m * Bt + ... , (4.2) 

and similar equations for the fermions. Here m 1B < m 2B < m 3B < m 4B . 

We now discuss how the contamination of higher excitations can be suppressed in (4.2), 
for the purpose of fitting the leading mode. In the limit of g — 0, we just have the simple 
harmonic oscillator with unit mass, frequency oo = vrt\ B , as well as fermionic excitations that 
are degenerate with the bosons. For g ^ 0, we expect a spectrum similar to the anharmonic 
oscillator, so that m 3B — m 1B > 2m 1B and similarly for the fermions. Thus the contribution of 
excited states to the Green function G 1B (t) can be neglected if t 2/mi B , since the relative 
exponential suppression in (4.2) will be something smaller than e~ 4 . To avoid finite size effects, 
which are of order exp[— (Na — t)m 1B ], we choose t ^ Na/3. (Recall that (3 = Na = 1 and 
mis > 10 in our simulations.) Taking these constraints on t into account, we can safely linearize 
the fit to obtain the leading excitation for each Green function (i.e., fit In G(t) to a linear function 
in t). 

Three different versions of the naive action have been simulated: S, eq. (2.8), which has 
no corrections; S c , eq. (3.8), which has the the 0(1) counterterm; and S ca , eq. (3.18), which 
includes the 0(a) improvement. 
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As a first case, we take g = 10, where the perturbation expansion parameter g/m 2 = 1/10 
is small. The results for m\B are shown in Fig. 3, and for m\p in Fig. 4. We compare these 
results to those of the continuum, obtained by a numerical solution of the Schrodinger equation 
corresponding to (2.1). It can be seen that the uncorrected naive action fails to approach the 
continuum limit for the boson. The discrepancy is of order 10 percent, consistent with an 
0(g/m 2 ) effect. The C(l)-subtracted and C(a)-improved results are in reasonable agreement 
with the continuum, though the latter is clearly better. The fermion mass does not appear to 
suffer from the same error, which is consistent with the fact that a counterterm is not required 
for it. However, it can be seen that the O (a) -improved action S ca more rapidly approaches the 
continuum result. (Recall that the correction to the fermion action occurs only at 0(a).) In 
fact, at the largest N that we were able to achieve, S ca obviously gives a superior approximation 
to the continuum results. We have also verified that the second set of excitations m 2 B and m^F 
approach the continuum results for the C(a)-improved action S ca . For brevity we do not present 
these results in detail. 




Figure 3: Leading boson mass for various forms of the naive action, with bare parameters m = g = 10. 
Large N corresponds to the continuum limit with (3 = Na held fixed at (5 = 1. Lines are drawn to 
guide the eye. 



Next, we have simulated the C(a)-improved action at strong coupling, g = 100. The 
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Figure 4: Similar to Fig. 3, except that the fermion mass is displayed. 



expansion parameter of perturbation theory would in this case be g/m 2 = 1, clearly beyond 
the range of validity of a perturbative approach. In Fig. 5 we give our results for the first 
boson and fermion excitations. It can be seen that, as indicated by the nonperturbative transfer 
matrix analysis of Section 3.2, the desired continuum limit is obtained. For comparison, we 
show the results of CG [2], obtained using the susy lattice action (2.18). Although the 0(a)- 
improvement of the naive lattice action gives a superior approximation to the continuum, an 
C(a)-improvement of the susy lattice action is straightforward to perform using the transfer 
matrix techniques of Section 3.2; this would also yield better agreement with the continuum. 



5. Conclusions 



In this article we have shown that a careful study of lattice power counting, or the transfer 
matrix, is necessary in order to understand the quantum continuum limit, even in the simple 
case of Id SQM. The lower UV degrees that occur in lattice perturbation theory with some 
amount of exact susy provide a compelling indication of the following: one should make exact 
lattice susy a requirement of the construction to the extent that it is possible. As has been seen, 
in the case of the susy lattice action of SQM, the continuum limit is obtained without the need 
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Figure 5: Leading boson and fermion mass for strong coupling, g = 100, m = 10, in the case of the 
O (a) -improved naive action. It can be seen that a good approximation to the continuum results is 
obtained for reasonable N. For comparison, we also plot the susy lattice results reported in [2] (the 
boson and fermion masses are degenerate). 



for counterterms, due to cancellations between diagrams with D = 0. The situation is similar in 
the field theory context (for example, compare to Appendix D of [3]); although, in field theory 
we expect that some amount of fine-tuning of lattice counterterms will generically be required, 
particularly if the continuum field theory is not finite. This is because modes near the UV cutoff 
will not decouple in the D > diagrams, yielding contributions that must be subtracted in 
order to match the continuum theory. Nevertheless, we expect the UV behavior of the lattice 
theory will be better if some exact susy is preserved. Practically speaking, this amounts to less 
counterterms that need to be calculated and adjusted. This is particularly important in a strong 
coupling regime, where the needed fine-tuning of lattice counterterms must be determined from 
a nonperturbative analysis such as Monte Carlo simulation. 

It is also of interest to study these issues in models where symmetries other than susy are 
partly responsible for good UV behavior in the continuum theory. For example, certain classes 
of nonlinear cr-models are known to be renormalizable in the continuum because the form of 
counterterms is greatly restricted by the nonlinearly realized symmetry. It has been found that it 



- 19 - 



is not possible to simultaneously preserve an exact lattice supersymmetry and the full nonlinear 
symmetry [3, 9]. Thus we expect the UV behavior of these models to be worse than in the 
continuum, and we are presently studying whether or not a finite set of counterterms suffices to 
guarantee the continuum limit. 
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Appendix 



A. The transfer matrix of the naive lattice partition function 



The purpose of this appendix is to show that eqn. (3.12) is the transfer matrix of the partition 
function of the naive lattice action (2.8), with Wilson parameter r = 1 (see (2.12)), supplemented 
for further use below by the counterterms (3.11): 



N 



S[k,£] = J2 



i=i 



la la 



, (A.l) 



where, as usual ki = k(xi),£i = £(xi), with i + N = i. The lattice partition function is then 
defined as: 



N 

Z[k,£] = cf[ I dft dp ! dx l e 

j=i „ 



(A.2) 



To continue, it is convenient to change the fermionic variables as follows: 

tpi = fji+i, fa = 7}i, 
so that the action (A.l) and partition function (A.2) become: 



(A.3) 



N 



^-(xi - Xi-if + ^-(h'ih'i + h) + rjiTji - r) i+1 r)i (1 + a(h" + £ t )) 
la la 



s[M = £ 

i=i 

Z[k,£] = cflj dffdrf J dx i e~ s , 



(A.4) 



where we absorbed a minus sign in the normalization constant c. 
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To construct the transfer matrix and hamiltonian (see, for example, [10]), we first introduce, 
at each time slice, a Hilbert space which is a tensor product of a bosonic and fermionic space. 
The bosonic Hilbert space is that of square integrable functions on the line. We use the basis 
of position eigenstates, {\x), (x'\x) = 5{x' — x)}, where the momentum and position operators, 
[p, q] = —i, act as q\x) = \x)x and e l ^ A \x) = \x + A) (note that we continue using the dimensions 
of the previous section: x has mass dimension —1/2, a has dimension of length, while the 
superpotential h(x) is dimensionless). The fermionic Hilbert space is two dimensional and is 
spanned by the vectors |0) and |1). The fermionic creation and annihilation operators obey 
{b',b} = 1, such that 6|0) = 0, |1) = 6*|0). The fermionic coherent states are defined as 
\ifj = |0) + |1)77, (v\ — (0| + where 77 and 77 are Grassmann variables. We then recall the 
usual relations for the decomposition of unity, (r)'\r]) = e^' n \ I = J dfjdrje^lr)) (r)\, and for traces 
of operators O on the fermionic Hilbert space: Tr O = J drfdrje'^ 11 (t)\0\ — 77); Tr(— 1) f O = 
f dfjdrje^ 11 (rj\0\rj) , with (— 1) F |0) = |0). We then define the transfer matrix by the equality: 

N 00 

Z[k,£] = c f[ I dffdrf I dx i e~ s = Tr (-1) F f N = (A.5) 
i=1 J -L 

N 00 

= f[ I dffdrf e'^ j dx l x 

x (r ) N ,x N \f\r ) N -\x N - 1 )(r ] N -\x N - 1 \f\r ) N - 2 ,x N - 2 ) x ... 
x(ri 2 ,x 2 \f \r ] 1 ,x 1 }(7 1 \x 1 \f\7 1 N ,x N ) , 

or, equivalently, through its matrix elements, (r]' t+1 , x l+1 \T\rj l , x l ) = 



~c l ' N exp 



i+l 



2 



.(A.6) 



(Substituting (A.6) into (A.5) one can immediately obtain (A. 4)). 

Now, we can use the fermion coherent state identity (following from the definitions given 
above) (77' |1 — XWb\rj) = e( 1 ~ x ^ r '' ri , as well as the action of p, q on the states \x l ) , to verify that 
the transfer matrix operator T, with matrix elements given by eqn. (A.6), is: 

00 

f = c^ N j dze^^-^ a - a -{h\q) 2 + k{q))^ exp (izp) (l + a(h"(q) + t(q)) Pb) ,(A.7) 
—00 

as already stated in eqn. (3.12) in the main text (where c 1 ^ is called M{a)). 

We note that the properties of the lattice partition function and transfer matrix of the 
super symmetric lattice action are studied in [3]. 
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